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<^ \ Abstract 

The exact two-particle energy eigenstates in a generic asymmetric rectangular box 
with periodic boundary conditions in all three directions are studied. Their relation 
with the elastic scattering phases of the two particles in the continuum are obtained 
. for both 1)4 and D2 symmetry. These results can be viewed as a generalization of 

\ the corresponding formulae in a cubic box obtained by Liischer before. In particular, 

the s-wave scattering length is related to the energy shift in the finite box. Possible 
. applications of these formulae are also discussed. 
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^ , 1 Introduction 



Scattering experiment serves as a major experimental tool in the study of 
interactions among particles. In these experiments, scattering cross sections 
are measured. By a partial wave analysis, one obtains the experimental re- 
sults on particle-particle scattering in terms scattering phase shifts in a chan- 
nel of definite quantum numbers. In the case of strong interaction, experi- 
mental results on hadron-hadron scattering phase shifts are available in the 
literature. On the theoretical side. Quantum Chromo dynamics (QCD) is be- 
lieved to be the underlying theory of strong interactions. However, due to 
its non-perturbative nature, low-energy hadron-hadron scattering should be 
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studied with a non-perturbative method. Lattice QCD provides a genuine 
non-perturbative method which can tackle these problems in principle, us- 
ing numerical simulations. In a typical lattice calculation, energy eigenvalues 
of two-particle states with definite symmetry can be obtained by measuring 
appropriate correlation functions. Therefore, it would be desirable to relate 
these energy eigenvalues which arc available through lattice calculations to the 
scattering phases which are obtained in the scattering experiment. This was 
accomplished in a series of papers by Liischer [1,2,3,4] for a cubic box topol- 
ogy. In these references, especially Ref. [3], Liischer found a non-perturbative 
relation of the energy of a two-particle state in a cubic box (a torus) with the 
corresponding elastic scattering phases of the two particles in the continuum. 
This formula, now known as Liischer's formula, has been utilized in a number 
of applications, e.g. linear sigma model in the broken phase [5], and also in 
quenched QCD [6,7,8,9,10,11,12,13]. Due to limited numerical computational 
power, the s-wave scattering length, which is related to the scattering phase 
shift at vanishing relative three momentum, is mostly studied in hadron scat- 
tering using quenched approximation. CP-PACS collaboration calculated the 
scattering phases at non-zero momenta in pion-pion s-wave scattering in the 
1 — 2 channel [12] using quenched Wilson fermions and recently also in two 
fiavor full QCD [14]. 

In typical lattice QCD calculations, if one would like to probe for physical in- 
formation concerning two-particle states with non-zero relative three momen- 
tum, large lattices have to be used which usually requires enormous amount of 
computing resource. One of the reasons for this difficulty is the following. In a 
cubic box, the three momenta of a single particle are quantized according to: 
k = (27r/L)n = (27r/L)(ni, n2, ^3), with n G Z^. ^ In order to control lattice 
artifacts due to these non-zero momentum modes, one needs to have large val- 
ues of L. One disadvantage of the cubic box is that the energy of a free particle 
with lowest non-zero momentum is degenerate. This means that the second 
lowest energy level of the particle with non- vanishing momentum corresponds 
to n = (1, 1, 0). If one would like to measure these states on the lattice, even 
larger values of L should be used. One way to remedy this is to use a three 
dimensional box whose shape is not cubic. If we use a generic rectangular box 
of size (rjiL) x {r)2L) x L with rji and 772 other than unity, we would have three 
different low-lying one-particle energy with non-zero momenta corresponding 
to n = (1,0,0), (0,1,0) and (0,0,1), respectively. This scenario is useful in 
practice since it presents more available low momentum modes for a given lat- 
tice size, which is important in the study of hadron-hadron scattering phase 
shift. Similar situation also occurs in the study of K to tttt matrix element (see 
Ref. [15] for a review and references therein). There, one also needs to study 
two-particle states with non- vanishing relative three-momentum. Again, a cu- 



^ We use the notation to stand for the set of three-dimensional integers. That 
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bic box yields too few available low-lying non-vanishing momenta and large 
value of L is needed to reach the physical interesting kinematic region. In all 
of these cases, one could try an asymmetric rectangular box with only one side 
being large while the other sides moderate. One only has to choose the pa- 
rameter 771 and 772 appropriately such that more low-lying momentum modes 
can be measure on the lattice with controllable lattice artifacts. 

In an asymmetric rectangular box, the original formulae due to Liischer, which 
give the relation between the energy eigenvalues of the two-particle states 
in the finite box and the continuum scattering phases, have to be modified 
accordingly. The purpose of this paper is to derive the equivalents of Liischer's 
formulae in the case of a generic rectangular (not necessarily cubic) box. 

We consider two-particle states in a box of size (r/iL) x (772 i^) x L with periodic 
boundary conditions in all three directions. For definiteness, we take r]i > 
1;^2 > 1, which amounts to denoting the the length of the smallest side of 
the rectangular box as L. The following derivation depends heavily on the 
previous results obtained in Ref. [3]. We will take over similar assumptions as 
in Ref. [3]. In particular, the relation between the energy eigenvalues and the 
scattering phases derived in the non-relativistic quantum mechanical model 
can be carried over to the case of relativistic, massive field theory under these 
assumptions, the same way as in the case of cubic box which was discussed 
in detail in Ref. [3]. For the quantum mechanical model, we assume that the 
range of the interaction, denoted by R, of the two-particle system is such that 
R < L/2. 

The modifications which have to be implemented, as compared with Ref. [3], 
are mainly concerned with different symmetries of the box. In a cubic box, 
the representations of the rotational group are decomposed into irreducible 
representations of the cubic group. In a generic asymmetric box, the symmetry 
of the system is reduced. In the case of r/i = 772 7^ 1, the basic group becomes 
1^4; if 771 7^ 772 7^ 1, the symmetry is further reduced to D2, modulo parity 
operation. Therefore, the final expression relating the energy eigenvalues of 
the system and the scattering phases will be different. 

This paper is organized as follows. In section 2, we discuss the singular periodic 
solutions to the Helmholtz equation. The energy eigenstates of the two-particle 
system can be expanded in terms of these solutions. In section 3, we discuss 
in detail the symmetry of an asymmetric box. Two cases are studied: 771 = 7/2 
in which case the basic symmetry group is D4 and 771 ^ 772 in which case 
the symmetry group is D2. The irreducible representations of the rotational 
group are decomposed into irreducible representations of these point groups. 
Energy eigenvalues in the A'l sector are related to the scattering phases for 
the two cases respectively. In section 4, we discuss the low-momentum and 
large-volume limit of the general formulae obtained in section 3. A simplified 
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formula is obtained for the scattering length and numerical values for the 
coefficients of this expansion arc listed. Finally in section 5, we conclude with 
some general remarks. Some details of the calculation are provided in the 
appendices for reference. 



Energy eigenstates and singular periodic solutions of Helmholtz 
equation 



Our notations close follow those used in Ref. [3]. The energy eigenstates in 
a periodic box is intimately related to the singular periodic solutions of the 
Helmholtz equation: 

(V^ + e)^lj(r) = . (1) 



If the function ip{r) is a solution to the Helmholtz equation for r 7^ 0; and it 
is periodic: '?/'(r + nL) = ipi'^)] ^-^d it is bounded by powers of near r = 0; 
we call '^(r) a singular periodic solution to the Helmholtz equation of degree 
A. 

The momentum modes in the rectangular box are qiiantized as: k = (27r/L)n. 
For every n = {rii, n2, ris) & Z^, we introduce the notations: 

n = (ni/771, 712/772, 713) , n = (7ii77i, 772712, 713) ■ (2) 



When discussing the singular periodic solutions to Helmholtz equation, one 
should differentiate two cases: regular values of k, which means that \k\ ^ 
(27r/L)|n| for any n G Z^' and singular values of k, which means that \k\ = 
(27r/L)|n| for some n G Z^'. For our purpose, it suffices to study the regular 
values of k. In this case, the singular periodic solutions of Helmholtz equation 
can be obtained from the Green's function: 



One can easily check that the function G{r; k^ is a singular periodic solution 
of Helmholtz equation with degree 1. More singular periodic solutions can be 
obtained as follows. We define: 

yi^r) = r%mm , (4) 



where Qr represents the solid angle parameters {6, (p) of r in spherical coor- 
dinates; Yi^ are the usual spherical harmonic functions. It is well-known that 
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3^im(r) consist of all linear independent, homogeneous functions in {x,y,z) of 
degree I that transform irreducibly under the rotational group. We then define: 

GUr;e)^yU^)G{r;e). (5) 



One can show that the functions Gimir; k"^) form a complete, linear indepen- 
dent set of functions of singular periodic solutions of the Helmholtz equation 
with degree /. That is to say, any singular periodic solution of the Helmholtz 
equation with degree A is given by 

V'W^E E vimGim{r,e), (6) 

1=0 m=-l 



with complex coefficients vim- The functions Gim{r; k"^) may be expanded into 
usual spherical harmonics with the result: 



Glm{v;k')^ 



2^ _ K-n^^' 



Air 



yim{^r)ni{kr) + J2 Mlm;l'm'Yl>m'{^r)jl'{kr) 
I'm' 



■(7) 



Here, ji and rii are the usual spherical Bcsscl functions and the matrix Aiim;i': 
is related to the modified zeta function via: 



\yl ^ {-yP ^Zl'm'il,q^;Vl,V2) //07 . lV07/^lV0o^1^ 

I'm' '/l'/2'' y 









^0 0^ 




m m' —s J 



X 



In this formula, the Wigner 3j-symbols can be related to the Clebcsh-Gordan 
coefficients in the usual way. [16] For a given angular momentum cutoff A, the 
quantity A4im;i'm' can be viewed as the matrix element of a linear operator M 
in a vector space Ti^, which is spanned by all harmonic polynomials of degree 
I < A. The modified zeta function is formally defined by: 

2im{s,q'^]Vi,V2) =Y.j^^^^^^-^ ■ (9) 



According to this definition, the modified zeta function at the right-hand side 
of Eq. (8) is formally divergent and needs to be analytically continued. Fol- 
lowing similar discussions as in Ref. [3], one could obtain a finite expression 
for the modified zeta function which is suitable for numerical evaluation. De- 
tailed formulae for Zim{s,q'^;r]i,r]2) at s = 1 and s — 2 are derived in the 
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appendix. Prom the analytically continued formula, it is obvious from the 
symmetry of D4 or D2 that, for / < 4, the only non-vanishing zeta functions 
at s = 1 arc: Zqq, Z20, 22±2, Z^q, Z/^±2 and Z/^±/^. It is also easy to verify that, 
if 771 = 772 = 1, all of the above definitions and formulae reduce to the those 
obtained in Ref . [3] . 

In the region where the interaction is vanishing, the energy eigenstates of the 
two-particle system can be expressed in terms of ordinary spherical Bessel 
functions: 

i}ira{r) = km [ai{k)ji{kr) + Pi{k)ni{kr)] , (10) 



for some constants bim- Also, this energy eigenfunction coincide with a singular 
periodic solution of the Helmholtz equation in this region. Comparison of 
Eq. (10) with Eq. (7) then yields: 



l'=0 m'=-l' 

blmm)-Vlm ^ \ ^ , (11) 

47r 

with Z = 0, 1, • • • , A. Note that this equation can be viewed as a linear equation 
in a vector space TYa, which is space of all complex vectors whose components 
are vim with / = 0, 1, ■ ■ ■ , A and m = —I, Matrix elements Aii'm';im can 

be viewed as the matrix element of an operator M in vector space TYa. The 
scattering phases of the two particles are related to the coefficients ai{k) and 
Pi{k) via: 



,2i5,(fc) _ aijk) + i(3i{k) 



ai{k) - i/3i{k) ' 



Therefore, non-trivial solution to Eq. (11) requires that: 



[12) 



det 



e^'' - U 



, {M + i)/{M -i) . (13) 



This gives the general relation between the energy eigenvalue of a two-particle 
eigenstate in a finite box with the corresponding scattering phases. 
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3 Symmetry of an asymmetric box 

The general result (13) obtained in the previous section can be further sim- 
plified when we consider irreducible representations of the symmetry group 
of the box. We know that energy eigenstates in a box can be characterized 



Table 1 

Basis polynomials in terms of spherical harmonics for various irreducible represen- 
tations r of the symmetry group 1)4 and D2 up to angular momentum I = 4. 





Group Di 


Group D2 


I 


r 


Basis polynomials 


r 


Basis polynomials 





At 


|0) 


A+ 




|0) = 


3^00 








1 


^2 


|1) 


Bi 




1 i \ 

|1) = 


3^10 










E~ 


(|i>, 1 1)) 


B2; B^ 


|i; 


— iJ'll + yi-i), 


1 \ 


[yii 


— yi- 


-1) 


2 


^1 


1 o\ 

|2) 






lo\ 

|2) = 


3^20 












\£}) 


R+. R+ 

B3 , B^ 


|2) 


= (3^21+3^2-1); 


2) = 


(3^21 


-3^2- 


-1) 




Bt; Bt 


|2);|2) 


A+; Bt 


|2) 


= (3^22+3^2-2); 


2) = 


(3^22 


-3^2- 


-2) 


3 




|3) 


BT 




|3) = 


3^30 










E- 


(|3),|3)) 


B2; B^ 


|3) 


= (3^31+3^3-1); 


3) = 


(3^31 


-3^3- 


-1) 




B2; Bi 


|3); 13) 


sr; A- 


|3) 


= (3^32 + 3^3-2); 


3) = 


{y32 


-3^3- 


-2) 




E- 


(|3),|3)) 


B2; B^ 


|3) 


= (3^33+3^3-3); 


3) = 


(3^33 


-3^3- 


-3) 


4 


At 


|4) 


A+ 




|4) = 


3^40 










E+ 


(|4),|4)) 


Bt; Bt 


|4) 


= (3^41 + 3^4-1); 


4) = 


(3^41 


-3^4- 


-1) 




Bt; Bt 


|4);|4) 


A+- Bt 


|4) 


= (3^42 + 3^4-2); 


4) = 


(3^42 


-3^4- 


-2) 




E+ 


(|4),|4)) 


Bt; Bt 


|4) 


= (3^43 + 3^4-3); 


4) = 


(3^43 


-3^4- 


-3) 




^1 ) ^2 


14); 1 4) 


A+; Bt 




|4) 


= (3^44 + 3^4-4); 


l)= 


(3^44 


-3^4- 


-4) 



by their transformation properties under the symmetry group of the box. For 
this purpose, one has to decompose the representations of the rotational group 
with angular momentum / into irreducible representations of the correspond- 
ing symmetry group of the box. For an asymmetric box, the relevant symmetry 
group is either D4 if 771 = 7^2 7^ 1, or if ?7i 7^ ?72 7^ 1- In a given symme- 
try sector, denoted by its irreducible representation F, the representation of 
the rotational group with angular momentum / is decomposed into irreducible 
representations of D4 or D2- This decomposition may contain the irreducible 
representation F. We may pick our basis as: |F, a; /, n). Here a runs from 1 to 
dim{T), the dimension of the irreducible representation F. Label n runs from 
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1 to the total mimber of occurrence of F in the decomposition of rotational 
group representation with angular momentum /. The matrix M is diagonal 
with respect to F and a by Schur's lemma. For convenience, we have listed 
these basis polynomials in terms of spherical harmonics yim{^) in Table 1. The 
corresponding shorthand notation for these basis are also listed. 

For the two-particle eigenstate in the symmetry sector F in a box of particular 
symmetry (either or D2), the energy eigenvalue, E = /c^/(2/i) with /j, being 
the reduced mass of the two particles, is determined by: 

det[e^'^ - U{r)] - , U{r) - (M(F) + i)/(M(F) - i) . (14) 

Here M(F) represents a linear operator in the vector space T-C\(T). This vector 
space is spanned by all complex vectors whose components are vin, with / < A, 
and n runs from 1 to the number of occurrence of F in the decomposition of 
representation with angular momentum I, see Ref. [3] for details. To write out 
more explicit formulae, one therefore has to consider decompositions of the 
rotational group representations under appropriate symmetries. 

We first consider the case rji = ri2- The basic symmetry group is D4, which 
has 4 one-dimensional (irreducible) representations: Ai, A2, Si, B2 and a 
two-dimensional irreducible representation E. ^The representations of the 
rotational group are decomposed according to: 

1 = + , 

2 = A+ + B+ + 3+ + E+ , (15) 

3 = ^2 + + +E- + E- , 

4 = A+ + A+ + A+ + B+ + B+ + E+ + E+ , 

In most lattice calculations, the symmetry sector that is easiest to investigate 
is the invariant sector: Af. We therefore will focus on this particular symmetry 
sector. As is seen, up to / < 4, s-wave, rf-wave and g wave contribute to this 
sector. This corresponds to /our linearly independent, homogeneous polynomi- 
als with degrees not more than 4, which are invariant under D4. Prom Table 1 
we see that the four polynomials can be identified as the basis: |0), |2), |4) and 

o 

|4), respectively. Therefore, we can write out the four-dimensional reduced 
matrix AI(A^) whose matrix elements are denoted as: Ai{Ai )ii' = mui — rrii'i, 

o 

with / and / takes values in 0, 2, 4 and 4, respectively. Using the general for- 
mula (8), it is straightforward to work out these reduced matrix elements in 

^ The notations of the irreducible representations of group D4 and D2 that we 
adopt here is taken from Ref. [17], chapter XII. 

^ Our convention for the spherical harmonics are taken from Ref. [18]. 
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terms of matrix elements A4im;i'm'- These are given explicitly in appendix B. 
We find that, in the case of D4 symmetry, Eq. (14) becomes: 



cot 5o — moo 


mo2 mo4 


171 
04 


m2o 


cot S2 — m22 m24 


m 
24 


m4o 


17142 cot 84 — 17144 


m 
44 


mo 

40 


mo mo 
42 44 


cot 84 — moo 

44 



= . (16) 



For the case rji 7^ r]2, the symmetry group becomes D2 which has only 4 
onc-dimcnsional irreducible representations: A, Bi, B2 and B^. The decompo- 
sition (15) is replaced by: 



= A+ , 

1^B^ + B2+B^ , 

2 = A+ + A+ + S+ + S+ + S3+ , (17) 

3= A- + + B^ + B2 + B2 + S3- + S3- , 

4 = + + + 5+ + S+ + Bt + + S+ + Bt , 



So, up to / < 4, occurs six times: once in / = and twice in / = 2 and 
three times in / = 4. The corresponding basis polynomials can be taken as: 

_ _ o ^ 

|0), |2), |2), |4), |4) and |4), respectively. The reduced matrix M{A^) is six- 
dimensional with matrix elements Aiw — mw. The explicit expressions for 
these matrix elements can by found in appendix B. The relation between the 
energy eigenvalue and the scattering phase is similar to Eq. (16) except that 
the matrix becomes a 6 x 6 matrix: 



- moo 


mo2 


mo2 


mo4 


moi 


m 

04 


m2o 


6 - ^22 


17122 


m24 


m24 


m 
24 


^20 


m22 


6 - ^22 


m24 


m24 


m_o 
24 


m4o 


m42 


m42 


^4 - m44 


m44 


m 
44 


mio 


m42 


m42 


m44 


^4 - ^44 


m_o 
44 


mo 

40 


mo 
42 


mo_ 
42 


mo 
44 


mo_ 
44 


^4 - moo 
44 



= , (18) 



where we have used the simplified notation: ^i{q) = cot 6i{q). In principle, if 
even higher angular momentum are desired, similar formulae can be derived. 
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4 Low momentum expansion and the scattering length 



Usually, scattering phases with higher angular momentum are much smaller 
than phases with lower angular momentum. This is particularly true in low- 
momentum scattering. It is well-known that for small relative momentum k, 
the scattering phases behave like: 

ta^n 5i(k) aik'^^+\ (19) 



for small k. Therefore, we anticipate that, in the low-momentum limit, scat- 
tering phases with small / will dominate the scattering process. If we treat the 
(i-wave and (7-wave scattering phases as small perturbations, we find that for 
the D4 symmetry, Eq. (16) can be simplified to: 

cot So - moo = + + . (20) 

cot 62 — m22 cot 04 — 17144 COt 64 — 17X44 



For the symmetry D2, similar to Eq. (20), Eq. (18) reads: 



2 2 
cot ^0 - moo = + 



cot 52 — m22 cot 82 — m22 

cot ^4 — m44 cot 64 — mil cot ^4 — mn 

If the (i-wave and g wave phase shift were small enough, it is easy to check 
that both Eq (20) and Eq (21) simphfies to: 

cot do(fc) = moo = ■ (22) 



For the general case, Eq. (20) and Eq. (21) offer the desired relation between 
the energy eigenvalues in the sector and the scattering phases for the cases 
rji = ri2 and ri\ 7^ 772, respectively. It is easy to verify that, in both Eq. (20) 
and Eq. (21), contributions that appear in the right-hand side of the equations 
are smaller by a factor of (f' compared with moo on the left-hand side. They 
are negligible as long as the relative momentum q is small enough. Therefore, 
in both cases, the s-wave scattering length will be determined by the zero 
momentum limit of Eq. (22). 

It is also possible to work out the corrections due to higher scattering phases 
to the s-wave scattering phase. For example, we have: 

n-K - 5o{q) = (j){q) + a2{q) tan 52(g) + (74(g) tan54(g) , (23) 
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Fig. 1. The functions a2{(^) and (74 (g^) as a function of are plotted for both D4 
symmetry with parameters ?7i = 772 = 2. 
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Fig. 2. The functions cr2(g'^) and iT4(g^) as a function of are plotted for both D2 
symmetry with parameters 771 = 2, 772 = 1-5. 

where the angle 0(g) is defined to via: — tan0(g) = l/moo(g)- The functions 
(72 (g) and £74(5) represents the sensitivity of higher scattering phases. For D4 
and D2 symmetry, they are given by: 
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■ "4/(1 +"^oo) . for group 

(J2(g)=<; (24) 

{ml2 + ml^) / {I + ml^) , for group D2. 



f^4(g) 



(mg4 + mg4)/(l + mlo) , for group D4. ^^^^ 

(m§4 + mg^ + ml~^)/{l + mgo) , for group D2. 



The functions cr2{q'^) and (74(5^) can be calculated using the matrix elements 
given in Appendix B. In Fig. 1 and Fig. 2, these functions are plotted versus 
for the case of D4 symmetry [rji = r/2 = 2) and for the case of D2 with 771 = 2.0, 
r]2 = 1.5, respectively. It is seen that the functions cr2(g^) and £74(5^) remain 
finite for all > 0. For some particular values of q^, however, these functions 
can become quite large in magnitude. This is due to almost coincidence of 
singularities of the numerator and denominator in matrix elements ruoi which 
happens for some choices of rji and 772- For values of away from these values, 
the functional values of (T2 and remain moderate. In Fig. 2, the lower two 
panels in the plot are simply the same function as in the upper panels with 
the scale of the vertical axis being magnified, in order to show the detailed 
variation of the functions. 

We remark here that, in principle, the corrections due to scattering phases 
with higher / can be estimated from lattice calculations as well. From Table 1 
it is seen that, for lattices with 1^4 symmetry, by inspecting energy eigenstates 
with Bi or B2 symmetry on the lattice, one can get an estimate for the 
d-wave scattering phase 82 which dominates these symmetry sectors. Similarly, 
for lattices with D2 symmetry, the eigenstates with Bi, B2 or B^ symme- 
try should be studied. It is also interesting to note that for lattices with D4 
symmetry, if we study the eigenstates with A2 symmetry, then the leading 
contribution comes from g^-wave scattering phase ^4. 

Wc now come to the discussion of scattering length. For a large box, a large 
L expansion of the formulae can be deduced. Using Eq. (22) and following 
similar derivations as in Ref. [3], we find that the s-wave scattering length ao 
is related to the energy difference in a generic rectangular box via: 



5E = 2^^° 



r]ir]2iJiL^ 



1 + ci{r]i, r]2) (^^) + 02(771, 772) ( 



L 



+ 



(26) 



Here, /i designates the reduced mass of the two particles whose mass values 
are tt^i and 1712, respectively. Energy shift 6E = E — rrii — 7712 where E is the 
energy eigenvalue of the two-particle state. Functions 01(771,7/2) and 02(771,772) 
are given by: 
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Zoo(l,0; 771,772) 

7r77i772 

^00(1, 0; Vi, V2) - ^00(2, 0; 771, 772) 

(7r77i772)2 



(27) 



where the subtracted zeta function is defined as: 



Zoo{s,q'';m,ri2)^ ^^2 _ ^2). • 



(28) 



Table 2 

Numerical values for the subtracted zeta functions and the coefHcients 01(771,172) 
and 02(7/1,772) under some typical topology. The three-dimensional rectangular box 

has a size Li = rjiL, L2 = 7?2^ and L3 = L. 



Li : L2 




m 


V2 


^00(1, 0; 771,7/2) 


■^00(2, 0; 7/1, 7/2) 


ci (7/1, 7/2) 


C2{m,V2) 


1:1: 


1 


1 


1 


-8.913633 


16.532316 


-2.837297 


6.375183 


6:5: 


4 


1.5 


1.25 


-12.964476 


41.526870 


-2.200918 


3.647224 


4:3: 


2 


2 


1.5 


-16.015122 


91.235227 


-1.699257 


1.860357 


3:2: 


2 


1.5 


1 


-10.974332 


32.259457 


-2.328826 


3.970732 


2:1: 


1 


2 


1 


-11.346631 


63.015304 


-1.805872 


1.664979 


3:3: 


2 


1.5 


1.5 


-14.430365 


53.784051 


-2.041479 


3.091200 


2:2: 


1 


2 


2 


-18.430516 


137.771800 


-1.466654 


1.278623 



In Table 2, we have listed numerical values for the coefficients 01(7/1,7/2) and 
C2{VijV2) under some typical topology. In the first column of the table, we 
tabulated the ratio for the three sides of the box: 771 : 772 : 1. Note that for 
^1 = ''72 = Ij these two functions reduce to the old numerical values for the 
cubic box which had been used in earlier scattering length calculations. 



It can be shown that contaminations from higher angular momentum scatter- 
ing phases come in at even higher powers of 1/L. For low relative momenta, 
the d-wave scattering phase behaves like: tan52{q) ~ 02^^ = ot2(27r/L)^g^, 
with a2 being the rf-wave scattering length. If we treat the effects due to tan52 
pcrturbatively, we see from Eq. (20) and Eq. (21) that, in Eq. (26), functions 
Ci and C2 receive contributions that are proportional to {a2/L^), which is of 
higher order in 1/L for large L. 
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5 Conclusions 



In this paper, we have studied two-particle scattering states in a generic rect- 
angular box with periodic boundary conditions. The relations of the energy 
eigenvalues and the scattering phases in the continuum arc found. These for- 
mulae can be viewed as a generalization of the well-known Liischer's formu- 
lae. In particular, we show that the s-wave scattering length is related to the 
energy shift by a simple formula, which is a direct generalization of the cor- 
responding formula in the case of cubic box. We argued that this asymmetric 
topology might be useful in practice since it provides more available low-lying 
momentum modes in a finite box, which will be advantageous in the study of 
scattering phase shifts at non-zero three momenta in hadron-hadron scattering 
and possibly also in other applications. 



A Appendix A 



In this appendix, some explicit formulae for the modified zeta function defined 
in Eq. (9) will be given. For convenience of analytic continuation, one first 
defines the heat kernel: 



Here the first equahty is the definition of the heat kernel /C(t, x) while the 
second follows from Poisson's identity. The summations which appear in this 
formula runs over all three-dimensional integers n E and the notations n 
and n are defined as in Eq. (2). It is evident that the heat kernel /C(t, x) is 
periodic in x with period 27rn, i.e. /C(t,x -|- 27rn) = /C(t,x). Given a positive 
number A > 0, we define the truncated heat kernel /C^(t, x) as: 



}C% x) = /C(t, x) - — ^ E e^'-^"*'' • (A.2) 

y"^^) |n|<A 



We may apply the operator 3^zm(— '^Vx) to the heat kernels defined above. 
These are denoted as: 

ICimit, X) = yim{-tV^)}C{t, X) , JCtit, X) = yim{-tV^)}C''{t, X) , (A.3) 
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Using heat kernels defined above, one could rewrite the modified zeta function 
as: 



Jdtf-'e"^']Ct{t,0). (A.4) 



This expression is convergent as long as Re{s) > (/ + 3)/2. Close to t = 0, the 
function e*^ ^tmi^^ 0) behaves like: 



(A.5) 



We may therefore analytically continue the zeta function by: 



(27r) 



'^l{)(\n{)lllll2 



2im{s,q';Vi,V2)-E^^^,_^,y ■ r(.) [{Anns -3/2) 

SloSmOlll1l2 



+ 



(47r) 2^3/2 



7 



(A.6) 



which is a vahd expression as long as Re{s) > 1/2. This completes the process 
of analytic continuation for the modified zeta functions. Using the explicit 
expression for K-i^it, 0), the integral from to 1 in the above formula can be 
further simplified. After some manipulations, we find that the results, when 
combined with the finite summation, is convergent for all A. Therefore, we 
can now send A to infinity and drop the third integral. Since we are only 
concerned with modified zeta functions at s = 1 or s = 2, we will only give 
the expressions for these two cases. The final expression may be written as: 



2^im{s,q^;Vi,V2) = e'''Yl 

n 

7^VlV2 



+ 



(2s - 3) 



(n^-q^y (n2_^2)s-i 



1 

SioSmo + ^SioSmoJ dtf-'>'\e'^' - 1) 



1 



n^O ^ 



(A.7) 



Note that this expression is vahd only for s = 1 or s = 2. 
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B Appendix B 



Table B.l 



Reduced matrix elements Mii'(T) for the symmetry group D4 for irreducible repre- 
sentations , A2 and E~. Here we have adopted the basis defined in Table 1. For 
all representations, only the non-vanishing matrix elements up to Z = 4 are listed. 



r 


|/) 


|/') 


yVl//,(r) 


At 








Woo 




2 





-W20 




2 


2 


Woo + ^W20 + f W40 




4 





W40 




4 


2 


-f W20 - ^W4o - ^T^Weo 




4 


4 


Woo + '"yf W20 + W40 + 'I^Weo + ''24^W80 






4 





^(W44 + W4-4) 






4 


2 


^(W44 + W4-4) - n715(W64 + W6-4) 






4 


4 


^(W44 + W4-4) - ^(W64 + W6-4) + ^^tHi^(W84 + W8-4) 






4 




4 


Woo 7l >^20 + 143 W40 u Jt3 + 143^/17^80 + (^§8 + ^^-^^ 


A — 

^2 


1 


1 


Woo + ;^W2o 




3 


1 


-^^20 - 75tW40 




3 


3 


Woo + 3j3W20+fiW40 + 33^J°3W60 


IT' — 


i 


i 


°° ~ 75*^20 T :^(W22 + W2-2j 




3 


3 


Woo + ;^W2o + ^W4o - HTj^Weo 

=F^(W22 + W2-2) T #(W42 + W4-2) T ^(W62 + W6-2) 




3 


i 


+ VI^^O ^ 2^(^22 + ^2-2) ^ 2^0^42 + W4-2) 




3 


3 


Woo - ^W2o + ^W4o - 55715 Weo T 7=5(W66 + We-e) 




3 


i 


-^(W22 + W2-2) + i73T(W42 + W4-2) ± 75(W44 + W4-4) 




3 


3 


-^(W22 + W2-2) + #(W42 + W4-2) - 55^(W62 + W6-2) 
±#(W44 + W4-4) T ^(W64 + W6-4) 
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Table B.2 

Reduced matrix elements A4u'(T) for the symmetry group D4 for irreducible repre- 
sentations E+, Bf, Bf, and A^. 



10 



\n 



E+ 



Woo + ^>V20 - |W40 T #(>V22 + W2-2) T 4^(W42 + W4-2) 

-40W20-^>V40 + ^W60 
^^(W22 + W2-2) ± ^(W42 + W4-2) ± (W62 + W6-2) 



Woo + + ^W40 



jWeo 



392 



77 "^^u ' 1001'"*" 11V13''"" 143V17 
T^(W22 + W2-2) T ^#(^42 + W4-2) 



Wso 



77 

/T05 



1001 

2iyi40, 



,,^,3(>Va2 + >V6_2)Tfi|^(>V82 + W. 



8-2) 



^(W22 + W2-2) - |^(W42 + W4-2) + ^(W62 + W6-2) 
±#(>V44 + W4^4) ± 1^(^64 + W6-4) 

Woo - ^>V2o - ^W4o + HTisWeo - j^mo 
T^(W66 + We-e) T ^(Wse + Wg-e) 

-^(W22 + W2-2) + J^,m2 + W4-2) 
+ ^(>^62 + W6-2) - tM?(W82 + W8-2) 



±^(W44 + W4-4) ± n^(W64 + W6-4) T Tf^(W84 + Ws 



13V187' 



B 



1/2 



Woo - ^W20 + yW40 ± ^(W44 + W4-4) 
Woo + ^W20 - i?W40 - 4?W60 + 



196 
143 v^7 



Wso 



11V13 



13V187 ' 



-4^W20 + ^W40 ,,^,3 



5^3 



Weo 



±^(W44 + W4-4) T ^(^64 + W6-4 



2/1 



Woo - ^W4o + iTTTjWeo ± ^(W44 + W4-4) ± ^(W64 + W6-4 



11V13 



At 



Woo - #W20 + ^W40 



TITE^^O + j£^mo - ^(W88 + W8-8) 



In this appendix, we list the reduced matrix elements that appeared in the 
formulae in the main text. We define: 



7^3/2^^^2g;+i 



(B.l) 



Using this notation, we now proceed to list relevant reduced matrix elements 
■Min-i'n'(r) for a given symmetry sector F. For the case of D4 symmetry, up to 
angular momentum / = 4, the non- vanishing matrix elements J^in;i'n'{^) are 
listed in Table B.l, Table B.2. 
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For the case of symmetry, by comparing with Table 1, one easily sees that 
the corresponding matrix elements are in fact the same as listed in Table B.l 
and Table B.2. The only difference is that the name of the irreducible repre- 
sentations are changed. For example, the matrix elements listed under E~ for 
symmetry are in fact corresponding matrix elements for the representation 
B2 and for the D2 symmetry, as suggested by Table 1. 
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